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Q\ \ Abstract 

£T) . The equilibrium of a hot dense plasma in a gravitational field is 

considered. From the standard equilibrium equations, the energy min- 
^ ■ imum at density about 10 25 particles per cm 3 and temperature about 

O '. 10 7 K was found. This effect plays an important role for astrophysics. 

It enables to explain the mechanism of the star magnetic field gener- 
' ation and to make a prediction for the spectrum of a star mass with 

• a wholly satisfactory agreement for the observation data. 



PACS: 64.30. +i: 95.30.-k 



1 The current hypothesis 



Now it is conventionally accepted to think that a density and temperature 
of a star interior substance are growing depthward of stars and can amount 
to tremendous values at their central cores. It seems that this growing is a 
necessary condition of an equilibrium of a star substance in the self gravity 
field. 

The substance exists as a hot dense electron-nuclear plasma at high pres- 
sure and temperature of a star interior. At this condition in zero approxima- 
tion, plasma can be considered as a Boltzmann ideal gas with energy 

E = -kTN (1) 
2 v ; 
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where N is the particle number. 

As the a direct inter-nuclear interaction is small in plasma, it can be 
neglected and one can write the equilibrium equation in the form [1]: 



/i e + m'lj) = const 



(2) 



where \i e is the electron chemical potential, ml = ^/n e is the mass of 
substance related to one electron, n e is the electron gas density, 7 is the 
mass density of plasma, ip is the Newton gravitational potential. Because 
Aip = —AnG'j, in a spherically symmetric case, it is reduced to 



is the function of particle density and temperature only. A single-meaning 
conclusion can be deduced from the equilibrium equation (Eq.(3)) if to con- 
sider hot non-degenerate plasma as an ideal gas: the balance of plasma par- 
ticles in the self gravity field demands a temperature and density increasing 
depthward of a star. 

The consideration of a hot dense electron-nuclear plasma as ideal gas is 
possible in zero approximation only because it means that the inter-particle 
interaction is completely neglected by definition. At first sight it is accept- 
able as the inter-particle interaction are small in comparison with the ideal 
gas energy (Eq.(l)). But the allowance of this interaction has a principal 
importance because it forms a stable equilibrium state of a hot plasma. It is 
lost at the consideration of plasma in the ideal gas approximation. 

2 The density and temperature of a hot dense 
plasmas in the equilibrium state 

2.1 The steady-state density of a hot dense plasma 

There are two main characteristic features which we must take to account at 
hot dense plasma consideration. The first of them is related to the quantum 




(3) 



The chemical potential of electron gas in the Boltzmann case [1] 
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properties of electron gas. The second feature is concerned with the presence 
of positively charged nuclei inside an electron gas of plasma. 



2.1.1 The correction for Fermi-statistics 

The estimation of electron gas energy in the Boltzmann case (kT ^> E F ) can 
be obtained by an expansion in series of the non-relativist ic Fermi-particle 
system full energy [1]: 

2 l/2 Vm 3/2 x £ 3/2 d 

E = 



J e (e-He)/kT + I W 



Where e is energy of particle. Because, in the Boltzmann case, \x e < 
and \fi e /kT\ 3> 1, the integrand at e^l kT <c 1 can be expanded into a series 
according to their powers e ^/ kT - £ / kT . 

As a result, the full energy of a hot electron gas [1] with account of its 
quantum properties is 



E = -kT 
2 



1 7r 3 / 2 /aoe 2 \ 3 / 2 



(6) 



where a = is Bohr radius. 

u me' 

It is important to underline that the correction for Fermi-statistic is posi- 
tive, because it takes into account that electron can not take places which are 
occupied by other electrons and a resulting pressure is more than pressure of 
an ideal gas at just the same density and temperature. 



2.1.2 The correction for a correlation of charged particles in plasma 

At high temperature, the plasma particles approach to an uniform space 
distribution. At this limit, the energy of ion and electron interaction is equal 
to zero. Some correlation in space distribution of particles arises as much as 
the positively charged particle groups itself preferably around particles with 
negative charges and vice versa. It is accepted to estimate the energy of this 
correlation by the method developed by Debye-Hwkkel for strong electrolytes 
[1]. The energy of a charged particle inside plasma is equal to e(p, where e is a 
charge of a particle, and ip is the electric potential induced by other particles 
on the particle under consideration. 

This potential inside plasma is determined by the Debye law [1]: 
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£ r_ 

(p(r) = -e r o (7) 



where the Debye radius is 



rD = y4^z (8) 



For small values of ratio the potential can be expanded into a series 



<p(r) = -- — + ... (9) 
r r D 

The following terms are converted into zero at r — > 0. The first term of this 
series is the potential of the considered particle. The second term is a poten- 
tial induced by other particles of plasma on the charge under consideration. 
Therefore, the correlation energy of plasma is [1] 



kT \dfj, e J n,t. 



(10) 



where Z is the nuclear charge and the density of nuclei is n n = 

Because the chemical potential of the Boltzmann ideal gas is determined 
by (Eq.4), at high temperature 

due = kT 
dn e n e 

and 

A V2 



5E = 



As an attraction between unlike charges which are placed nearer one to an- 
other is prevalent over a repulsion of like charges, the plasma pressure is 
below the pressure of the ideal gas at the same parameters. For this reason 
this correction is a negative one. 
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2.1.3 The density of a hot dense plasma at the equilibrium state 

Thereby the full energy of plasma with account of the main corrections on 
inter-particle interaction is 



E = lkTN e 
2 



(13) 



At a constant full number of particles in the system and at constant 
temperature, the condition of the equilibrium state existence at the energy 
minimum is 

,an e J n,t 

what allows one to obtain the steady-state value of density of hot non- 
relativistic plasma 

n* = 16( Q Z 2 + 3 1)3 * 2 • 10 24 (Z + I) 3 cm-\ (15) 
97r 2 ao 

Here the Fermi-energy of electron gas of equilibrium plasma at this density 
is 

16\ 2//3 me 4 / „ . „ N 9 e 2 

a 



e '^= (ff¥< z+i ^ L 4 (z+ 1)2 (16) 



2.2 Equilibrium temperature of a hot non-relativist ic 
star 

As the steady-state value of density of hot non-relativistic plasma is known, 
we can obtain a steady-state value of temperature of a hot non-relativistic 
plasma. 

According to the virial theorem [1, 3], the potential energy U of particles 
with Coulomb interaction is equal to their double kinetic energy T with an 
opposite sign 

U = -2T (17) 
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and their full energy is equal to kinetic energy with an opposite sign. At 
neglecting of small corrections at high temperature, one can write the full 
energy of hot dense plasma as 

Eplasma = U + ^kTN = ~kTN. (18) 

As the plasma temperature is high enough, pressure of black radiation cannot 
be neglected. The full energy of a star depending on the plasma energy and 
the black radiation energy is 

3 7r 2 / kT\ 3 

^ = __*™ + _(_)vH" (19) 

The equilibrium temperature of a body consisting of hot non-relativistic 
plasma is determined by the energy minimum condition 

/ 9E tota i \ 
V dT ) N y 

It gives 

MSf (Z + 1) ^ K2 - (Z + 1) ' 107 *' (21) 

Usually all substances have a positive thermal capacity. Therefore the 
minimal energy for such substances exists at T = 0. The existence, in our 
case, of the energy minimum at finite temperature T* 7^ is not confusing. 
Each small part of a star has a positive thermal capacity, but a gravitational 
interaction of these parts between themselves results into that the thermal 
capacity of a star as the whole becomes negative at some temperature and 
a star energy decreases with an increased temperature. As a result there 
are two branches of the temperature dependence of a star energy - with 
the negative capacity at low temperature and with a positive capacity at 
high temperature. Between them, at some finite temperature T* there is a 
minimum of energy. 

Steady-state values of density and temperature for hot non-relativistic 
plasma have been considered above. One can see, that the criterion according 
to which plasma may be considered as a hot one 

kT > e F (22) 
6 



(20) 



with account of Eq.16 is satisfied if the nuclear charge is not too big 



kT* / 45 y/ 3 he 0.3 _! 40 

7^7) (z + i)e^YTi a w zTT' (23J 

where ct = e 2 //ic = 1/137 is the fine structure constant. 



3 The equilibrium of a dense electron-nuclear 
plasma 

According to definition (Eq.(4)) in plasma equilibrium state at a constant 
temperature and density, its chemical potential must be constant too: 

/i(n*,T*) = const. (24) 

(One can note, that for a cold degenerate plasma the chemical potential 
does not depend on temperature and it is constant at n — const.) 

Thus the equilibrium of plasma in gravity field can be obtained if the role 
of the field of another nature than gravitational one, for example, the electric 
field, is taken into account. Speaking more precisely the equilibrium equation 
(Eq.(2)) must imply all fields, which can have an impact on particles, for 
example, electric field for system of charged particles: 

f^ + Yl = const (25) 

i 

(where qi and fa are charge and potential of any nature (gravitational,electric)). 
In a spherical symmetric case at /i — const it reduces to 

G 7 m' = pq. (26) 
where p = q ■ n e is the electric charge density, 

q = G 1/2 m' (27) 

is the charge induced in plasma cell by the gravity (related to one electron). 
One should not think that the gravity field really induces in plasma some 
additional charge. We can rather speak about electric polarization of plasma 
which can be described as some redistribution of internal charges in the 
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plasma body. Essentially it stays electrically neutral as a whole, because the 
positive charge with volume density 



p = G 



1/2 7 

Z 



(28) 



is concentrated inside the charged plasma core and the corresponding nega- 
tive electric charge exists on its surface. 

Thus, in a spherical symmetric case the equilibrium equation can be 
rewritten as 



The Thomas-Fermi approximation permits to consider the balance in 
plasma cells in more detail [4]. 

3.0.1 The equilibrium density of another kind of dense plasmas 

The above consideration of equilibrium of hot non-relativistic plasma is ob- 
vious, but it is a characteristic property not of this plasma only. 

The direct consideration of a plasma equilibrium in Fermi-Thomas ap- 
proximation shows that the application of a gravity field to plasma induces 
its electric polarization ([4], [11]) and so the equilibrium equation in the 
form Eq.(29) is applicable to all kind of dense plasmas - relativistic or non- 
relativistic and simultaneously degenerate or non-degenerate ones. 

Essentially at that one can find the constancy of density and chemical 
potential of another kind of plasmas in an equilibrium condition. 

So for a cold non-relativistic plasma the kinetic energy of electron is 



7g + pE = 



(29) 



where 



E = 



g 



(30) 



(Gzy/ 2 ' 




(31) 



Its potential energy is 





(32) 
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As according to the virial theorem Ek ~ — E v the equilibrium electron density 
is 



n e ^a 3 . (33) 

and it does not depend on temperature. 

The relativistic plasma exists at huge pressure which induces the neu- 
tronization of substance. For this process density 

Ue = 37W (34) 
is characteristic [1]. Where A is the difference of nuclear bonding energy of 
neighbouring interacting nuclei. Because of this difference 

A « m e c 2 (35) 

the equilibrium of relativistic plasma density (at condition of a homogeneous 
mixing of reacting substance) 

n e ~ -i-(cm r 3 ~ 10 30 cm" 3 (36) 
on 2 

Where a = e 2 /he is the fine structure constant, a = h 2 /m e e 2 is the Bohr 
radius. 

A similar consideration can be extended on the neutron matter if it is 
considered as electron-proton plasma in the neutron environment. 



4 The giro- magnetic ratio of stars 

Gravitation produces a redistribution of free charges in plasma inside a star. 
Essentially, a star as a whole conserves its electric neutrality. However, as 
the star rotates about its axis, positive volume charges are moving on smaller 
radial distances than the surface negative charge. It induces a magnetic field 
which can be measured. 

The magnetic moment of surface spherical layer, which carries the charge 
Q, is 

/i_ = ~Q{IR 2 (37) 
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Where f2 is rotational velocity. 

The magnetic moment induced by a volume charge is 



^ + = ^li / pr'dr = l -QSlR 2 (38) 
J Jo 5 

Where Q = f pi? 3 . 

Thus, the summary magnetic moment of a star is 

^ = -^(yP^ 3 )0R 3 (39) 

As the mass density inside a star is constant, the angular momentum of 
a star is 

L = -MVLR 2 (40) 
5 

According to Eq.(27), we obtain that the giro-magnetic ratio is expressed 
through world constants only: 

* = £ = f (41 > 

It can be verified by measurement data. 

The values of giro-magnetic ratio for all celestial bodies (for which they 
are known today) are shown in Fig.l. 

The data for planets are taken from [6], the data for stars are taken from 
[7], and those for pulsars - from [8]. Therefore, for all celestial bodies - for 
planets and their satellites, for Ap-stars and several pulsars - the calculated 
value of gyromagnetic ratio Eq.(41) with a logarithmic accuracy quite sat- 
isfactorily agrees with measurements, which angular momenta themselves 
change within the limits of more than 20 orders. 



5 The stellar mass distribution 

5.1 The mass of star consisting of hot non-relativist ic 
dense plasma 

Inside the stellar core consisting of hot dense plasma the gravitational force 
is counterbalanced by the electric force (Eq.(29)) and gradient of pressure is 
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Figure 1: The observed values of the magnetic moments of celestial bodies 
vs. their angular momenta. On the ordinate, the logarithm of the magnetic 
moment over Gs ■ cm 3 is plotted; on the abscissa the logarithm of the angular 
momentum over erg ■ s is shown. The solid line illustrates Eq.(41). The dash- 
dotted line is the fitting of the observed values. 
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absent. The absence of a pressure gradient inside the hot star core does not 
mean absence of pressure. It is not difficult to make sure, if the gravity force 
inside a star is compensated by electric force, the negative energy of gravita- 
tional field is cancelled by the energy of electric field. The non-compensated 
part of energy is the energy of gravitational field outside a star. This field 
has an energy 

where M* and R* are the mass and radius of a hot plasma core. This external 
gravitational field tends to compress a star. 

By definition, pressure of a gas (at entropy S=0) is 

9E G \ _ _Eg = GMl_ 
dV ) s= o 3K 6iW 1 ' 

where V* = ^-R+ is the volume of a star. 

The equilibrium of a star exists when the gravitational pressure, which 
tends to compress a star, is counterbalanced by internal pressure of the hot 
electron gas and by pressure of the radiation: 

GMl 7T 2 (kTM 

W = lT «"' + 45W <44> 
Finally we have the mass of a star core consisting from hot dense plasma 

M - 1 ^M(^-X' 2 (-\ 2 m « 6A7Mch « 1^ (45) 

where Men = ( q^j J fn P = 3.42 • 10 33 g is the Chandrasechar mass. It is 

important to underline that the obtained stellar mass estimation Eq.(45) is 
depending on only one parameter AjZ. 

One can note that the equilibrium radius of star core 



_ (3/2)3 /loy/*/ hc y /2 ao 

R *-—r\-) \Grt) (z+i)a/z- (46) 



v 

such as depends on A and Z only. 
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5.2 The mass of star consisting of cold relativistic plasma 

With an increase of density, the plasma can turn into relativistic state. It 
occurs when Fermi momentum of electrons 

p F = (37T 2 ) 1/3 nl /3 fi > m e c (47) 

This value of momentum corresponds to the steady-state density of a 
substance under neutronization (Eq.(36)) at n* « 10 30 cm~ 3 . At temperature 

T < ^! » 10 10 K. (48) 

it can be considered as cold. 

As the pressure of a relativistic electron gas is 

Pr = ^f^n^hc, (49) 
in accordance with Eq.(43), the pressure balance obtains the form: 

£K = ( 50) 

Therefore, the relativistic degenerate star in equilibrium state must have 
a steady value of mass 

M _ j he \ 3 / 2 m p ^ 4.9M Cfe ^ 9.OM 

at radius corresponding to Eq.(36) 

(he y/ 2 aa ^ 10-*R Q 

The objects which have such masses and density are best suited to the 
dwarfs. 
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Figure 2: The mass distribution of binary stars [9]. On abscissa, the loga- 
rithm of the star mass over the Sun mass is shown. Solid lines mark masses 
which agree with selected values of A/Z from Eq.(45) for stars. The dotted 
lines mark masses which agree with selected values of A/Z from Eq.(51) for 
dwarfs. 



5.3 The comparison of calculated star masses with ob- 
servations 

The comparison of calculated results with the data of measurements is shown 
in Fig.2. There is an extensive star mass evidence but the mass measurement 
of binary stars has a sufficient accuracy only. The mass distribution of visual 
and eclipsing binary stars [9] is shown in Fig.2. On abscissa, the logarithm 
of the star mass over the Sun mass is plotted. Solid lines mark masses which 
agree with selected values of A/Z for stars from Eq.(45). The dotted lines 
mark A/Z for dwarfs from Eq.(51). 
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According to existing knowledge, the hydrogen inside dwarfs is fully burnt 
out. In full agreement with it in Fig. 2 there are hydrogen stars and there are 
not dwarfs with A/Z = 1, whereas there are both - helium-deuterium stars 
and dwarfs - with A/Z = 2. Attention should be attract to the fact that 
there is the peak of stars consisting of heavy nuclei with A/Z = 2.8 and peak 
for dwarfs with A/Z = 3 where the Sun is placed. 

The nuclei with A/Z > 2.8 are absent in terrestrial condition. Evidently 
a situation for stars composed of nuclei with A/Z > 3 is more complicated 
and it demands a special and more attentive consideration. 



5.4 The mass of a star consisting of neutron matter 

Dwarfs may be considered as stars where a process of neutronisation is just 
beginning. Finally, at nuclear density, plasma turns into neutron matter. 

It is agreed that a pulsar is a star consisting of neutron matter with some 
impurity of other particles - electrons and protons. Evidently, at nuclear 
density neutrons and protons are indistinguishable inside pulsars as inside a 
huge nucleus. The neutron matter can be considered as some kind of plasma 
where electrons and protons exist in a neutron environment. Under this 
condition a very small impurity of electrons and protons (about a level of 
1CT 18 ) is enough to induce a sufficient electric polarization and to balance 
the action of gravity. 

It should be noted that neutron matter at nuclear density is a degenerate 
substance, if T 10 12 K. It is important that the neutron gas is not ultra- 
relativistic at a nuclear density. It is a relativistic gas only when the Fermi 
momentum of neutrons is pp ~ m p c. For this reason, the equation of pressure 
balance for neutron relativistic gas has a complicated form [1]: 

G K^-^L( l - ahi -ls h Ui) (53) 



6RV 327r 2 ft 3 V3 s 3 2 
Where f = 4Arsh^. 

As a result, the mass of a pulsar is 



M pulsar = ( - J M Ch F(rf) (54) 



where 
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and 




ft(37r 2 7)V3 




(56) 



depends on world constants and the density of neutron matter 7 only. 

The estimation shows that at pp ~ m p c the function ^(7) 1. The 
numerical calculation of ^(7) shows that the calculated value of mass of a 
pulsar (Eq.(54)) fully coincides with observed data (Fig. (3)) when the matter 
density 7 approaches the nuclear density 3 • 10 14 g/cm? in full agreement 
with the assumption of incompressibility of nuclear matter. (At this density, 
« 0.4 and Ffa) « 0.72). 

The mass distribution of pulsars is shown in Fig. (3). In the upper scale 
the density of neutron matter according to Eqs.(54)-(56) is plotted. 
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0.3 0.4 
log M/Mo 



Figure 3: The mass distribution of pulsars [10]. In lower abscissa, the loga- 
rithm of the pulsar mass over the Sun mass is shown. In upper abscissa, the 
density of the substance in g ■ cm" 3 according to Eqs.(54)-(56) is plotted. 
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6 Conclusion. 



To conclude, it is important to note that no disputable speculative assump- 
tions were made above. All above results was obtained on the standard 
physical base by standard formal methods. 

The novelty of the developed approach to plasma equilibrium consists in 
the rejection of the usually accepted point of view that an increase of tem- 
perature and density depthward of a celestial body is a requirement of an 
equilibrium of a substance in a gravity field. It is really applicable to equi- 
librium of atomic substances. It is shown above that there is the minimum 
for electron-nuclear plasma energy at some density and temperature which 
brings to the equilibrium in a gravitation field at zero gradient of the general 
parameters of plasma because of its electrical polarization. 

The understanding of this effect discovers the simple mechanism of gen- 
eration of the magnetic field by celestial bodies. It can be noted that all 
previous models tried to solve the other basic problem: they tried to cal- 
culate the magnetic field of a celestial body. Now space flights and the 
development of astronomy discovered a remarkable and previously unknown 
fact: the magnetic moments of all celestial bodies are proportional to their 
angular momenta and the proportionality coefficient is determined by the 
ratio of world constants only ([5], [6]). The explanation of this phenomenon 
is really the basic problem of planetary and stars magnetism nowadays. The 
developed theory gives a simple and standard solution to this problem. 

This approach gives a possibility to predict important properties of stars 
in their steady state. Starting from the equilibrium conditions it allows 
to calculate masses of different types of stars. Thus the masses of stars 
composed by non-relativistic non-degenerate plasma and dwarfs composed 
by relativistic degenerate plasma can be expressed by the ratio of world 
constants and one variable parameter (A/Z) only, and this statement is in a 
rather good agreement with the observation data. Just as the predicted value 
of mass of pulsars is in full agreement with observations at the assumption 
of incompressibility of nuclear matter. 

Some considered questions are described more systematically in [4], [11]. 
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